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A theory of dynamical control by modulation for optimal decoherence reduction is developed. 
It is based on the non-Markovian Euler-Lagrange equation for the energy-constrained field that 
minimizes the average dephasing rate of a qubit for any given dephasing spectrum. 
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Introduction. — Quantum information processing 
(QIP) harbors enormous unleashed potential in the form 
of efficient algorithms for classically intractable tasks 
and unconditionally secure cryptography 1]. Perhaps 
the largest hurdle on the way to a realization of this 
potential is the problem of decoherence, which results 
when a quantum system, such as a quantum computer, 
interacts with an uncontrollable environment or bath 
[2|. Decoherence reduces the information processing 
capabilities of quantum computers to the point where 
they can be efficiently simulated on a classical computer. 
In spite of dramatic progress in the form of a theory 
of fault tolerant quantum error correction (QEC) [3], 
finding methods for overcoming decoherence that are 
both efficient and practical remains an important 
challenge. An alternative to QEC that is substantially 
less resource-intensive is dynamical decoupling (DD) 
0) S ul- In DD one applies a succession of short and 
strong pulses to the system, designed to stroboscopically 
decouple it from the environment. This can significantly 
slow down decoherence, though not halt it completely, 
since unlike QEC, DD does not contain an entropy 
removal mechanism. Similar in spirit to DD (in the 
sense of being feedback- free) , but more general, is the 
method we term here "dynamical control by modula- 
tion" (DCM), wherein one may apply to the system a 
sequence of arbitrarily-shaped pulses whose duration 
may vary anywhere from the stroboscopic limit to that of 
continuous dynamical modulation 7, 8, % [l(J, [ll|, LL2| • In 
the DCM approach, the decoherence rate is governed by 
a universal expression, in the form of an overlap between 
the bath-response and modulation spectra, subject to 
finite spectral bandwidth and amplitude constraints. 

Neither DD 0, i, [| nor DCM @, El studies have 
so far gone beyond particular schemes for suppression of 
decoherence. What is lacking is a systematic theory for 
finding the optimal modulation for any given decoher- 
ence process. Here we apply variational principles to the 
DCM approach in order to address this problem. We de- 
rive an equation for the optimal, energy- constrained con- 
trol by modulation (ODCM) that minimizes dephasing, 
for any given dephasing spectrum. We numerically solve 



this equation, and compare the optimal modulation to 
energy-constrained DD pulses. We show that ODCM al- 
ways outperforms DD when subjected to the same energy 
constraint. We note that Ref. Q developed an optimal 
DD pulse sequence for the diagonal spin-boson model of 
pure dephasing, but without an energy constraint, i.e., 
assuming zero-width pulses. This was improved upon by 
perturbatively accounting for pulse widths in Ref. [13j |. 

Model. — We consider a driven two-level system (qubit) 
with ground and excited states \g) and |e) separated by 
energy uj a (we set K= 1), and Hamiltonian 



H(t) = (u> a + 6 r (t))\e)(e\ + (V(t)\g)(e\ 



(1) 



where V(t) = Vt(i)e~ iuJat + c.c. is a time-dependent res- 
onant classical driving field with amplitude f2(£), and 
S r (u!) describes random, Gaussian distributed, zero-mean 
energy fluctuations. Let IV'(i)) denote the solution of the 
time-dependent Schrodinger equation with the Hamilto- 
nian H(t), and let the density matrix p(t) = \ip(t))(ip(t)\ 
denote the corresponding ensemble average over realiza- 
tions of 5 r (t). We are interested in the average fidelity 
(F(t)), where (• • • ) is the average over all possible initial 
pure states of the fidelity F(t) = \(yb(0)\p(t)\yb(0)}\. It 
can be shown that 1121: 



(F(t)) = 1 - aR(t)t 



(2) 



R(t) = 2Rc 



dfa*(ti-t 2 )c*(ti)c(t2)) (3) 



$(£) = S r (t)8 r (0) e(t) = e- i ti dt ^ tl) 



(4) 



where < a < 1 is a known constant, (•)£* = j J* -dti is 
the time-average, R(t) is the average modified dephasing 
rate, $(i) is the second ensemble-average moment of the 
random (stationary non-Markov) noise, and e(t) is the 
phase factor associated with the modulation. 
We impose the energy bound constraint 



Jo 



dt\n(t)\' 



E 



(5) 



where T is the total modulation time and E is the energy 
constraint. As a boundary condition we require that the 
field is turned on, i.e. O(0) = 0. 



Although the analysis below is given in the time- 
domain, it is advantageous to analyze the problem in the 
frequency domain, in terms of the universal expressions 

as 



R(t) = 2tt / dujG(uj)F t (uj) 



dt$(t)e l 



G(lu) = (2tt)- 



F t (w) = \e t (v)\ 2 /t e t (w) 



(6) 

(7) 



where G(u>) is the dephasing spectrum, et(u>) is the finite- 
time Fourier transform of the modulation function, and 
F t (u) is the normalized spectral modulation intensity. 

The model we have just described applies to a 
qubit undergoing dephasing due to coupling to a finite- 
temperature bath of harmonic oscillators with energies 
huj\. The qubit then has an average modified dephasing 
rate of the form given by Eqs. (O,© where the dephasing 
spectrum is given by [12j : 

G(w) = (n(w) + l)Go(w)+n(-w)G (-w) (9) 



Go(w) = V| K a| 2 <5(w-co a ) 



A 



(10) 



where Gq{lS) is the zero-temperature bath spectrum, n\ 
is the off-diagonal coupling coefficient of the qubit to the 
bath oscillator A, and n(ui) 



o/3w 



1 



is the aver- 



age number of quanta in the oscillator (bath mode) with 
frequency w, with /3 the inverse temperature. 

Since R(t) is the overlap between the dephasing and 
modulation spectra, it can be reduced by choosing an ap- 
propriate modulation that reduces this overlap [3, M, H2] • 
We shall show that the optimal modulation reduces the 
spectral overlap of the dephasing and modulation spec- 
tra (Fig. [5]). However, since the energy constraint in the 
frequency domain is non-trivial we shall derive the equa- 
tions for optimal modulation using the time domain. 

Optimization. — We wish to find the optimal modula- 
tion, i.e., time-dependent near-resonant field, that mini- 
mizes R(t). Calculus of variation is an often- used tech- 
nique of optimal control theory, e.g., [IjjLLjI. We apply 
it to derive the Euler-Lagrange (EL) equations for the 
energy-constrained optimal modulation. The accumu- 
lated phase due to the modulation is <f>(t) = L drCl(T). 

Let us write $(i) = $(i)e lAt , where $(£) and A are the 
amplitude and spectral center of the correlation func- 
tion, respectively. Using Eqs. @ and (O, we can then 
derive the EL equation for the optimal modulation (see 
Appendix [ft] and |B) : 

\4>(t) = -Z[t,<j>(t)] (11) 

Z[t,4>(t)} = (^{\t-ti\)wa[<f>(t)-<f>(ti)+A(t-ti)] 



where A is the Lagrange multiplier. The boundary con- 
ditions for the accumulated phase are 0(0) = 0(0) = 0, 



which results in a smooth solution and accounts for turn- 
ing the control field on. Eliminating A we find that the 
optimal control field shape is the solution to the following 
equation (see Appendix [Bj: 



m 



-VEZ[t,4>{t)\ 



\j ft dh\fc dt 2 Z[t 2 ,<)>{t 2 )] 



(12) 



Equation (|12p is the central result of this work. It fur- 
nishes the optimal time-dependent field amplitude, that 
maximizes the average fidelity (F(t)) at the final time 
T, via Q(t) = <p(t). Although Eq. (TT2"1) is a complicated 
non-linear integro-differential equation, it is very useful 
indeed, as we show next. 

Linearized EL equation. — Assuming that we have a 
good initial guess </>o(£) for the modulation, we can look 
for the optimal deviation v(t) by writing <p(t) — cf>o(t) + 
v(t), where v{t) <C 1. Expanding Eq. (jTTj) in powers of 
v(t) and retaining only the first order, the linearized EL 
equation becomes (see Appendix [Cjl : 

Xi)(t) + (Q(t, t i; 4> Q (t)) ("(*) ~ K*i)))t = -C(*i <h, A) 
Q(t,ti;&(t)) = *(|t-ti|)x 

cos(0o(t)-^o(ti) + A(t-ti)) 
C(t; 0o(t), A) = XMt) + Z[t, Mt)}- (13) 

This linearized EL equation is valid also in the case of 
short time optimal modulation, for which we simply set 
<t> (t) = 0, subject to v{t) < 1 for < t < T < 1. 

Numerical analysis. — Armed with the equations for 
the optimal modulation, we turn to solving them numer- 
ically for specific decoherence scenarios, defined by their 
dephasing spectra G{ui). We obtain the numerical solu- 
tion to the integro-differential Eq. (fT2|) via an iterative 
process, where we guess a probable solution that satisfies 
the boundary conditions and the constraint, use it in the 
RHS of Eq. (fT2|) to compute the integral, and solve the 
resulting differential equation. The solution is then used 
in the RHS of Eq. [[12")). and so on. 

For the examples presented below, we checked that sev- 
eral initial guesses converged to the same optimal mod- 
ulation. Most importantly, we found that the optimal 
modulation is robust against random control field imper- 
fections. This is due to the fact that the decoherence rate 
is determined by the accumulated phase and not the in- 
stantaneous modulation, Eq. ([J|. Specifically, we found 
that a 10% zero-mean random pulse fluctuation results 
in less than a 1% increase in the optimal dephasing rate. 

We compare the optimal dephasing rate to the one ob- 
tained by the popular periodic DD control ( "bang bang" ) 
procedure (J|, but to make the comparison meaningful 
we impose the same energy constraint. Finite-duration 
periodic DD against pure dephasing is the "bang bang" 
application of n 7r-pulses and is given in our setting by 

m) = l 7t/v jr - t<JT+ " j = °--- n - 1 (u) 

y ' I otherwise v ' 



where v < r is the width of each pulse and r is the in- 
terval between pulses. The energy constraint E and the 
total modulation duration T = tit + v are related via 
n = vEjii 2 . In the frequency domain, the spectral mod- 
ulation intensity can be described by a series of peaks, 
where the two main peaks are at ±7t/t. Thus, the peaks 
are shifted in proportion to the energy invested in the 
modulation. However, DD is not an admissible solution 
to our EL equation due to its discontinuous derivative. In 
order to improve the comparison, we apply our linearized 
EL equation with the DD modulation as an initial guess, 
and obtain the optimal modulation in the vicinity of the 
DD control. 

(a) Single-peak resonant dephasing spectrum,.- This 
simple dephasing spectrum describes a common scenario 
where $(£) = e~'/* c 7/t c , where 7 is the long-time de- 
phasing rate [R(t — > 00) = 2irj] and t c is the noise cor- 
relation time. Fig. [lja) shows R(T), normalized to the 
bare (unmodulated) dephasing rate, as a function of the 
energy constraint. As expected, the more energy is avail- 
able for modulation, the lower is the dephasing rate. For 
low energies the optimal modulation significantly outper- 
forms DD, while at higher energies this difference disap- 
pears. These results can be understood from Fig.[2]Ja), by 
noticing that the two central DD peaks have significant 
overlap with G(uS) at the low energy value shown. As E 
is increased at fixed T the DD peaks move farther apart, 
and have less overlap with G(u>), leading to improved 
performance. Applying the linearized EL equation with 
DD as initial guess yields only mild improvements (not 
shown). The explanation for the superior performance 
of the optimal modulation is also evident from Fig. (Ha): 
since higher frequencies have lower coupling strength in 
this case, the optimal control "reshapes" so as to maxi- 
mize its weight in the high-frequency range, to the extent 
permitted by the energy constraint. The modulation can 
be well approximated by Q(t) = o[l + e _t / T (t/T - 1)], 
where a is determined by the energy constraint, which 
fits the inset in Fig.QJa). 

(b) Single-peak off-resonant dephasing spectrum. — 
This dephasing spectrum describes a variation on the 
aforementioned scenario, where the spectral peak is 
shifted [A ^ in $(*) = <l(£)e iAt ], e.g., coupling to 
a non-resonant bath. With no other constraints, the 
optimal modulation is trivially similar to the one of 
the resonant spectrum, with a shifted energy-constraint 
^non-res _ j^rcs _|_ ^ However, by imposing a positivity 
constraint, <j>(t) > (positive field amplitude), one ob- 
tains non-trivial behavior of R(T) as a function of the 
energy constraint - see Fig. HJb). Here we used the lin- 
earized EL equation with the DD modulation (JT4J) as an 
initial guess. For both the DD and optimal modulations, 
we observe an initial increase in the dephasing rate as 
a function of energy, followed by a decrease. For DD, 
this can be interpreted as a manifestation of the initial 
anti-Zeno effect and the subsequent quantum Zeno effect 
[3, [lCJ- Because of the positivity constraint, the optimal 
modulation does worse than the unmodulated case, for 




FIG. 1: (color online) Average modified final decoherence rate 
R(T), normalized with respect to the unmodulated rate, as a 
function of energy constraint. DD - dash, magenta. Optimal 
modulation - solid, blue. Insets: optimal modulation Q(t) for 
different energy constraints, (a) Single-peak resonant dephas- 
ing spectrum (inset: E = 20). (b) Single-peak off-resonant 
spectrum (inset: E = 50). (c) 1// spectrum (inset: E = 30). 
(d) Multi-peaked spectrum (inset: E = 30). 
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FIG. 2: (color online) Dephasing spectrum G(uj) (solid, red), 
optimal (dot-dash, green) and DD (dash, magenta) modula- 
tion spectra Ft{oj), in arbitrary units (a.u.). Same parame- 
ters as in the insets of Fig. [TJ 



low enough energy. The DD modulation is optimal for 
small energy constraints, hence the decoherence rates of 
DD and our optimal solution coincide. This is because 
the DD peaks do not overlap the off-resonant spectral 
peak. However, as the positive-frequency main DD peak 
[Fig. Elb)] nears the off-resonant spectral peak, with in- 
creased energy, the optimal modulation diverges from the 
DD modulation, and "reshapes" itself so as to couple to 
higher modes of the bath. In the time domain [Fig. [TJb) 
inset] , this is seen as a smoothing of the abrupt DD mod- 
ulation. At even higher energy constraints, there is once 
more no improvement by the optimal modulation over 
DD, yet there is an improvement over the unmodulated 
case. Over the entire range of E, the optimal modula- 



tion results in a much natter R(T) than DD, which is an 
indicator of its robustness. While DD is strongly influ- 
enced by the off-resonant peak, the optimal modulation 
exploits the energy available to find the minimal overlap, 
irrespective of dephasing spectrum. 

(c) 1// dephasing spectrum. — The ubiquitous 1// de- 
phasing spectrum that describes a variety of experiments 
- e.g., charge noise in superconducting qubits [16| - is 
given in our notation by G(u>) oc 1/w, with cutoffs u> m i n 
and w max . Fig.[T]Jc) shows that as expected, the more en- 
ergy is available for modulation, the lower the dephasing 
rate. Since, as in case (a), higher frequencies now have 
lower coupling strength, the optimal control "reshapes" 
so as to have as high a weight in the high frequency 
range as the energy constraint allows [Fig. Etc)]. This 
is expressed in the time-domain [Fig. [TJc) inset] as the 
initial increase in the modulation strength (t < 50). The 
later decrease in modulation strength can be attributed 
to the lower cutoff, where the optimal modulation bene- 
fits from lower frequencies, i.e., lower modulation ampli- 
tudes. Upon comparing the 1// case to the Lorentzian 
spectrum, Fig. [TJa) , we observe a similar optimal initial 
chirped modulation in the time domain. Despite the dif- 
ferences in the long-time behavior (due to the lower cutoff 
in the 1// case), these two examples allow us to gener- 
alize to any dephasing spectrum with a monotonically 
decreasing system-bath coupling strength as a function 
of frequency. The optimal modulation for such spectra 
will be an energy-constrained chirped modulation, with 
variations due to other spectral characteristics, e.g., cut- 
offs. 

(d) Multi-peaked dephasing spectrum. — This describes 
the most general scenario, where there can be several 



resonances and noise correlation times. Fig. [TJd) shows 
R(T) as a function of the energy constraint. Once again, 
because DD does not account for the dephasing spec- 
trum, its performance is much worse than the optimal 
modulation, whose "reshaping" results in monotonically 
improving performance: the peaks of the optimal modu- 
lation are predominantly anti-correlated with the peaks 
of G{uj). 

Conclusions. — We have found the optimal modulation 
for countering pure dephasing upon imposing an energy 
constraint on the DCM approach 0, [l2| , by deriving and 
solving the Euler-Lagrange equation p^]) . This yields op- 
timal reduction of the overlap of the dephasing and the 
modulation intensity spectra. We stress that our opti- 
mal control theory results are also applicable to scenarios 
other than pure dephasing, such as amplitude noise (re- 
laxation), due to the universality of Eqs. ©-([H]) [71 KL2J]. 
The form of the energy constraint will then differ in de- 
tail from the pure dephasing case. However, our general 
conclusions about the optimal modulation to minimize 
spectral overlap, will remain valid. We expect that the 
optimal modulation technique will find useful applica- 
tions in quantum information processing and quantum 
computation. The price is that one must acquire inti- 
mate knowledge of the noise spectrum, which is often ne- 
glected, as previous control techniques such as DD and 
QEC had no use for it. We have shown that this infor- 
mation can result in the maximization of fidelity, under 
operational constraints. 
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APPENDIX A: GENERAL COMMENTS ON DERIVING OPTIMAL FUNCTIONS 



For the optimal control of a functional 



with the constraint 



one follows the following procedure: 
(i) Solve the Euler-Lagrange equation: 



F(y,y) = f dtF{t,y,y) 

Jo 

K{y,y) = I dtK(t,y,y) = E 

Jo 



(Al) 



(A2) 



SF d SF _ 
Sy dt Sy 



5K d SK 
Sy dt Sy 



(A3) 



where SF is the variation of F, A is the Lagrange multiplier, and the boundary conditions are y(0) = yo and j)(0) = y\. 
(ii) Insert the solution y(i; A) into the constraint: 



lC(y(t;\),y(t;\))=E 



(A4) 



and obtain A = X(E). 

(iii) Eliminate A by inserting X(E) into y(t; A) and obtain the optimal solution, y(t; E), that minimizes the functional 
F, under the constraint K, = E. 



APPENDIX B: DERIVATION OF THE EULER-LAGRANGE EQUATION 

The average modified decoherence rate is given by: 



R(T) 



dt dti®(t - h) cos[(j>(t) - (f)(ti) + A(t - ti)] 
T Jo Jo 



2 
T 



T rT 

dt / dt x Q{t - ti)$(t - h) cos[0(i) - <t>{t x ) + A{t - h)] 

x Jo 



where &(t) is the Heaviside step function. 

One arrives at the following variation of the average modified decoherence rate: 



5R{T) = 1 f dt [ dt x 
1 Jo Jo 

-Q(t - ti)$(t - ii) sin[0(i) - <j>{t\) + A(t - h)]S<j)(t) 
+Q{t - ti)$(t - h) sm[<j)(t) - 4>{ti) + A(t - ti)]^(*i) 



, dt / dh 
F Jo Jo 

-Q(t - ti)$(t - ti) sm[<p(t) - cj)(ti) + A(t - h)}6<j>(t) 
-0(ti - t)&(h - t) $m[cj>(t) - <f>(ti) + A(t - h)]S<f>(t) 

-$(|t - ii|) sin^(i) - 4>{ti) + A(t - h)]6(f>(t) 



= - I dt J dh 



o Jo 



(Bl) 
(B2) 



(B3) 



(B4) 
(B5) 



where we have made a t <-> t\ substitution in the second integrand of Eq. (|B3[) . and notice that @(t)f(t)+Q(—t)f(—t) 
f(\t\). 



One can easily see that denning the constraint functional as: 



K(t,<j>{t),<Kt))= / dh\4>{t 1 )\ 2 = E, 
Jo 

with K(t,4>(t),(j)(t)) = \(f>(ti)\ 2 results in the variation: 

5K = 2<t>(t)8<j)(t). 

Combining Eqs. (|B5|) . (|B7|) and (|A3|) results in the Euler-Lagrange equation: 

\4>(t) + Z[t,cf>(t)] =0 

where 

Z[t, (f>(t)] = - / dti*(|* - ti|) sh#(i) - ^(ti) + A(i - ti)]. 



(B6) 

(B7) 
(B8) 

(B9) 



APPENDIX C: DERIVATION OF THE LINEARIZED EULER-LAGRANGE EQUATIONS 



In some cases it is advantageous to linearize the EL equations with respect to the modulation. If one looks for 
the optimal deviation v{t) from a given pulse shape, 4>o(t), then one can write <p{t) = 4>o(t) + ^(t), where u(t) <C 1. 
Equation (|B8[) then becomes: 



X(Mt) + v(t)) 



1 
TJ 



dti«(|t - ti|) sin[0 o (t) - o (ti) + A(t - ti)] coa[v(t) - i/(*i)] 



+ cos[^o(i) - Mh) + A(t - ti)] sin[z/(*) - v{h)] 



dti«(|t-ti|) sin[^ o (t)-0o(ti) + A(t-ti)] 
J Jo L 

+ cos[0 o (t) - 0o (*i) + A(t - ti)][i/(t) - i/(ti)]l + 0(i/ 2 (t)) 



(CI) 



(C2) 



where we approximated sin[i/(t) — ^(ti)] s=s ^(t) — ^(ti) + 0(i/ 3 (t)) and cos[^(t) — f(<i)] ~ 1 — j(^(t) — ^(ii)) 2 + 0(i/ 4 (i)). 
The linearized Euler-Lagrange equation becomes: 



where 



**>(*) + ™ / *iO(*, *i; &(*)) ("(*) - K*0) = -c(*; 0o> A ) 

J Jo 



Q(t,h;Mt)) = *(|*-*i|)oo8(^o(t)-0o(*i) + A(t-ti)) 

C(t;0 o (t),A) = A0o(*)+Z[t,0o(t)]. 

I 



(C3) 



(C4) 
(C5) 



